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ABSTRACT 

We discuss the submillimeter properties of Extremely Red Objects (EROs) in the two Canada-UK Deep Sub- 
millimeter Survey (CUDSS) fields. We measure the mean submillimeter flux of the ERO population (to K'< 
20.7) and find 0.40 ± 0.07 mJy for EROs selected by (I-K) > 4.0 and 0.56 ± 0.09 mJy for EROs selected by 
(R — K) > 5.3 but, these measurements are dominated by discrete, bright submillimeter sources. We estimate 
that EROs produ ce 7- 1 1 % of the far-infrare d background at 850/im. This is substantially less than a previous 
measurement bv lWehner. Barger. & Kneib (2002) and we discuss possible reasons for this discrepancy. We 
show that ERO counterparts to bright submillimeter sources lie within the starburst region of the near-infrared 
color-color plot of Pozzetti & Mannuccij (120001) . Finally, we claim that pairs or small groups of EROs with 
separations of < 10 arcsec often mark regions of strong submillimeter flux. 

Subject headings: galaxies:evolution-galaxies:formation-galaxies:high-redshift-galaxies:starburst- 
infrared:galaxies-submillimeter 



1. INTRODUCTION 

A complete accounting of all the systems which produce 
the extragalactic background light is key to our understanding 
of galaxy assembly and evolution. Observations of the early 
universe have revealed a wide variety of objects whose rela- 
tionship to each other, importance to galaxy formation, and 
contribution to the extragalactic background energy remains 
unclear. The fact that these systems have been selected over 
different wavelength ranges and using different selection cri- 
teria has made it is difficult to determine if we are studying 
different aspects of the same objects, or completely orthogo- 
nal populations. Many of these objects are found over simi- 
lar redshift ranges and are co-evolving, while others may be 
linked through time by evolutionary processes. 

At near-infrared wavelengths a population of red ob- 
jects begins to emerge at f a int, K > 18, magnitudes 
(lElston. Rieke. & Riekd Il988t IMcCarthv. Persson. & Westl 
1992; Hu & Ridgwavl ll994l) . This class of systems, labeled 
Extremely Red Objects (EROs), is photometrically defined by 
NIR color, typically (I-K) > 4 or (R-K) > 5 - 5.3. Two dif- 
ferent galaxy populations are selected by this color definition: 
an old stellar population at z ~ 1 , or alternatively, a young star- 
burst at high-redshift, deeply enshrouded in dust. Low-mass 
stars and brown-dwarfs within the Galaxy also have similar 
NIR colors but have such low surface densities relative to faint 
galaxies that they do not significantly contaminate ERO sam- 
ples. 

Determining the fraction EROs belonging to each category 
has been difficult work, hindered by the faintness of these 
objects at optical and NIR wavelengths. Spectroscopic ob- 
servations indicate that evolved stellar populations and dusty 

1 Leiden Observatory, Niels Bohrweg 2, 2333 CA Leiden, The Nether- 
lands 

2 Visiting Astronomer, Canada-France-Hawaii Telescope, Operated by the 
National Research Council of Canada, the Centre de la Recherche Scien- 
tifique de France, and the University of Hawaii 

4 Department of Astronomy and Astrophysics, University of Toronto, 60 
St George St, Toronto, Ontario, Canada, M5S 1A1 

3 Department of Physics and Astronomy, Cardiff University, RO. Box 9 1 3. 
Cardiff, CF2 3YB, UK 

5 Department of Physics, Swiss Federal Institute of Technology (ETH- 
Zurich), ETH Hoenggerberg, CH-8093, Zurich, Switzerland 



starbu rsts are present in roughly equal num bers dDunk^p^etalJ 
1 19961 [ Graham & Devil 1 9961 ICohen et aljl!999t Omatti et alJ 
2002a) and this picture is generally supported by larger stud- 
ies employing less ex pensive, but more ambiguo u s, tech- 
niques such as color fPozzetti & Mannucc ij |2QQQ} |Mar|:inil 
1200 lt lS mail et al. 2002a), radio emission ( Smail et al. 2002a), 
or morphology ( Yan & Thompson 2003; Gilbank et al. 2004). 
There is substantial variation between the results how- 
ever, stemming from different survey dep ths, varying color- 
selection criteri a, and clustering e ffects (|Daddi et al J | 2000t 
I McCarthv et al] 1200 It Isiriith et alJl2002t I Smail et al.l l2002a: 
Dadd i et all2002tlYan & Thompsonl2003l) . 

The detection of submillimeter emission from an ERO 
is a direct confirmation of a starburst nature (ignoring 
for the moment the problem of confusion within the 
beam). At current observational limits however submillime- 
ter imaging can only detect the most luminous star form- 
ing EROs, those that overlap with the discrete submillimeter- 
bright population at hig h-redshi f t discovered with S CUBA 
(SMGs) J^ail. Iyison. & Blainl fl997t iBargeretalJ [19981 
Hugh es et alJl998tlEales et alJ2000tlScott et alJ2002l) which 
are sites of perhaps the most substantial star formation in 
the high-redshift universe. Though still suffering from small 
number statistics ^20-30% of the systems detected deep sub- 



millimeter surveys are EROs 


when observed in the NIR) 


( Smail et alj|1999|: llvison et alJ 


2002: ChaDmaTetaT] r 2003a: 



Webb et al. 2003c) though this fraction obviously depends 
on the relative depths of the the 850/.tm and /f-band imag- 
ing and the uncertainty in the in the identifications of the 
counterparts of non-radio detected SMGs. What is not clear, 
however, is how many of the ^-selected EROs are discrete, 
bright submillimeter sources and the submillimeter propeties 
of EROs below the current detection limits of submillimeter 
surveys. These less luminous systems, though not individu- 
ally submillimeter-bright could together produce a substan- 
tial fraction of the FIR background light. Two previous stud- 
ies have attemp ted to address these issues through statistical 
studie s EROs ( iMohan et al1l2002t IWehner. Barger. & Kneibl 
2002) and the results of both groups indicate that as a pop- 
ulation the se systems emit substantia l submillimeter flux. In 
particular, Wehner. Barger. & Kneit] d2002h (hereafter W02) 
found K' < 2 1 .6 EROs produce ~50% of the FIR background 
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light at 850/im, with this contribution peaking at K' ^20. 
However, as with all ERO studies, these results may depend 
on ERO selection criteria and field-to-field variance. 

Thus, the relationship between EROs and bright SMGs, 
and the submillimeter properties of EROs below the de- 
tection limit of SCUBA, remain open questions, as does 
the ERO contribution to the far-infrared background energy. 
Those EROs which are extremely luminous in the submil- 
limeter are a minority of both the ERO population and the 
SMGs (at least to current K depths) but must represent an 
important phase of galaxy evolution. Understanding their 
natures will provide a key link between these two popula- 
tions. Isolating the properties of EROs which are individu- 
ally or statistically submillimeter bright which differentiate 
them from the rest of the ERO population is an important 
step toward this understanding. Such properties may include 
the NIR colors and magnitudes ( Pozzetti & Mannucci l2000t 
IWehner. Barger. & K neib 2002]), morphologies, redshifts and 
the local environment. Environment plays a key role in 
galaxy evolution and there is growing evidence of a corre- 
lation between submillimeter emission and high-dens ity re- 
gions or galaxy merging and intera ctions Jlvison et alJl2000t 
Smail et al.ll2003l lAlexander et al.ll2003l iWebb et alJl2003 d 
Neri et alJ2003l] Chapman et al. 2003c) and this may be be an 
important driver of star-formation in the ERO population. Ad- 
ditionally, the different redshift distributions of the EROs and 
SMGs (as far as we know them) (Chapman et al. 2003bJ hint 
at an evolutionary relationship between these populations per- 
haps with SMGs fading into less luminous star-forming EROs 
at later times. Addressing these issues will not only illuminate 
the natures of these two populations but the process of galaxy 
evolution in general. 

The paper is organized as follows. In §2 and §3 we present 
the data and define our ERO samples. In §4 we measure the 
measure submillimeter flux of EROs. We discuss the ERO 
contribution to the ex tragalacti c background light in §5. In §6 
we present the Pozzetti & Mannucci (2000) color-color plot 
for these systems. In §7 we discuss the effect of ERO compan- 
ions in the analysis. We finish in §8 with some thoughts differ- 
ence between submillimeter-bright and faint EROs. Through- 
out we have used a flat, A=0.7 cosmology and H = 72 
km/s/Mpc. 

2. OBSERVATIONS 

2.1. Submillimeter Data 

The Canada-UK Deep Submillimeter Survey (CUDSS) was 
a deep, random-field submillimeter survey carried out from 
1998-2002 at 850^m using SCUBA on the JCMT The sub- 
millimeter data a nd catalogs are discuss ed extensively in 
lEales et alJ feOOOl) and Webb et al. ( 2003b). The survey con- 
sists of two fields, CUDSS 14 and CUDSS03, which cover 48 
arcmin 2 and 60 arcmin 2 respectively, to a roughly uniform 3a 
depth of 3 mJy. The chopping and nodding technique em- 
ployed by SCUBA to remove sky-flux creates distinct beam 
pattern on the maps: each positive source is flanked by two 
negative "sources" with half the flux-density, offset from the 
positive source by a distance equal to the chop throw (in our 
survey this is 30")- This results in a map which has a total flux 
of zero. In addition to the bright source catalog (i.e. S$so^ m > 
3 mJy) the contiguous nature of t he maps allow an e xtraction 
information to fainter flux levels (Web betall2003al) . 

2.2. Optical and NIR Imaging 



The two CUDSS fields were originally selected to lie within 
Canada-France Redshift Survey fields dLillv etalJ fl995) and 
are well studied at many wavelengths. The Canada-France 
Deep Fields Survey has obtained n ew deep optical imag- 
ing, including RIz, o ver these areas JMcCraken et al.ll200ll 
Brod win et all 12004) using the CFHT 12K camera. The 
pixel scale of these data is 0.207 arcsec/pixel, and reaches 
2a depths of 7=25.58 (25.06), R=25.12 (25.72), and z=24.55 
(24.75) for CUDSS03 (CUDSS 14). 

The K' data were obtained in 2000 with the CFHT-IR 
camera on the Canad a-France-Hawaii Telescope (CFHT) 
JHutchings et all2002tlWebb et all2003d) . The data cover ~ 
2/3 of the two CUDSS fields, and areas of 35.4' and 55.8' for 
CUDSS 14 and CUDSS03 respectively. The CFHT-IR camera 
field-of-view is 3.6' x 3.6' and the final images are mosaics of 
many individual pointings. The pixel scale is identical to that 
of the RIz data, 0.207 arcseconds/pixel, and conditions dur- 
ing the run were stable with seeing typically 0.8". The depth 
over the two fields is not uniform but a 5a limit of K'= 20.7 
is reached over most of the area. 

3. ERO SAMPLE SELECTION 

Object detection was performed on the K' images usin g the 
source extraction package SExtractor ( B ertin & Arnoutl 19961) 
and magnitudes were extracted through 3" aperture photom- 
etry. The K' source counts are presented in Fig^and are in 
good agreement with other surveys, though for clarity we plot 
only the K20 counts (Cimatti et al. 2002b). Colors were mea- 
sured in 2.5" apertures on images degraded to common see- 
ing of 1.3". The effects o f galactic extinction were removed 
in both fields JMcCraken et al.l200ll) 

We selected EROs by two criteria commonly used in the 
literature, (I-K) > 4 and (R-K) > 5.3. For objects without 
R or I detections or detections at < 2a (as determined from 
the average noise properties of the R and I images) the 2a 
limit on their magnitude was adopted. This provides a lower 
(blue) limit on their color, and they were taken to be EROs 
if this limit met the above selection criteria. We present the 
ERO source counts (as a function of K' magnitude) in Fig. [2] 
and compare them with other surveys which apply the same 
selection criteria. We find the source counts of each field, 
using both color criteria, are in good agreement with these 
other surveys and therefore conclude that these fields do not 
contain exceptional over or under-densities of EROs. 

In CUDSS03 (CUDSS 14) we select 130 (1 1 1) objects with 
(I-K) > 4 and 77 (87) with (R-K) > 5.3. There are 51 and 77 
galaxies in common between the two samples for CUDSS03 
and CUDSS 14 respectively. Thus, (I-K) > 4.0 appears to be 
a more generous ERO definition than (R- K) > 5.3, an effect 
also seen bv lSmail et alJ J2002al) . However, the reddest col- 
ors we can measure for our faintest K magnitudes are (I-K) = 
5.0 ((I-K) = 4.6) for CUDSS03 (CUDSS 14) and (R-K) = 5.2 
for both fields and therefore we probe EROs to slightly deeper 
K-magnitudes using (7-/T)-selection. 

4. THE SUBMILLIMETER FLUX OF EROS 

4.1. The overlap of the ERO and discrete submillimeter 
populations 

A number of recent papers have investigated the near- 
infrared properties of SMGs which have been r o bustly 
identi fied through radio detections (llvison et alJ 2002; 
Chap man et al.ll2003at IWebb et al.H20Q3d) . While these sys- 
tems show a range of near-infrared colors, and some are com- 
posite systems with components of different colors, the av- 
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FIG. 1.— The K number counts for both fields, CUDSS 14 and 
CUDSS 03. Shown for com parison are the number counts from the K20 Sur- 
vey. ICimatti et aTl2 002b ) 
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FIG. 2. — The number counts for EROs selected by (I-K) > 4.0 and (R- 
K) > 5.3 for both fields, C UDSS 14 and CUDS S03. Shown fo r comparison 
are the number counts from Roche. Dunlop, & Almaini 1 2003), Smith et al. 
1 2002), i Daddi et al. 2000), and Smail et al. 1 2002jJT(see legend). 



erage (I-K) color for these systems is significantly redder 
than that of the background ^-selected population, and many 
SMGs are classified as EROs. 

In the CUDSS catalog there are ten bright (discrete) 
SMGs (5850/im > 3mJy) whose most likely counterpart is 
an ERO (T able [TJ and th ese objects are discuss ed in detail 
in lWebb et all d2003d) and lClements et al.1 (120031) . These pa- 
pers also describe the statistical identification methods which 
were used to select these counterparts, but we briefly outline 
them again here. We worked with radio (1.4 GHz) and IK 
imaging and selected the statistically most likely counterparts 
based on positional probability analysis using surface densi- 
ties determined from eith er the 1 .4 GHz data or / and K mag- 
nitudes and colors. In Webb et al. (2003c) we showed that 
the IK identification algorithm blindly and unambiguously se- 
lected the same counterpart as the radio data. We concluded 
that if the SMG source properties are similar to the radio- 



identified systems at redshifts beyond the reach of the radio, 
the IK identification algorithm will remain successful at se- 
lecting the correct counterparts. Three of the ten EROs are 
detected at 1 .4GHz which allows us to state with further con- 
fidence that these particular identifications are correct. The 
remaining seven, however, are not detected in the radio (per- 
haps because they lie at redshifts beyond the radio limit) and 
so, although statistically very likely, some (or all) of these 
identifications could be incorrect. Thus, between 6-20% of 
the discrete SMGs in the CUDSS catalogue have ERO coun- 
terparts to K = 20.7. 

4.2. The mean submillimeter flux of the EROs 

To estimate the mean 850/im flux of the ERO population we 
measure the flux on the SC UBA map at t he A"'-s elected posi- 
tion of each ERO (see also Webb et al. (2003a)). The com- 
bination of the large beam size of the JCMT at 850/im (15" 
FWHM) and the chopping technique (30" in EW) means the 
measured flux of an ERO could be contaminated by a nearby 
bright submillimeter source which is unrelated to the ERO. 
For example, if the location of an ERO fell within the negative 
off-position of a bright submillimeter source we would mea- 
sure an erroneous negative flux for the ERO. To reduce this ef- 
fect we have removed all of the discrete submillimeter sources 
from the maps which are securely identifi ed (i.e. ide ntified 
through radio-detecti ons) with non- EROs (lEales et alJ feOQO: 
Webb et al.ll2003blcl IClements et"alll2003l) before extracting 
the ERO fluxes. We then co mpute the noise-w ei ghted mean 
of these measurements (see lEales et al] (120001) : IWebb et all 
(2003b} for noise discussion). These results are listed in 
Table |2] In practice, removing the discrete submillimeter 
sources which are securely identified with non-EROs has not 
significantly (i.e. within the uncertainties) changed the mean 
flux measurements. It has lead to an increase in the mea- 
sured mean flux of 0.12mJy for CUDSS03 and a decrease of 
0.05mJy for CUDSS 14 (remembering the change can be ei- 
ther negative or positive because of the negative-positive na- 
ture of the beam profile) and thus the combined measurement 
has remained constant within the uncertainties. 

Though we use the true -determined position for the 
EROs, the submillimeter-determined position (if it were pos- 
sible to measure) will a lways be offset from this location du e 
to noise and confusion (Eales et al. 2000; Webb et al. 2003c). 
At the ^'-position there is an equal chance of noise boosting 
the signal up or down, while if we attempt to search for the 
peak flux of th ERO we would be pulled toward positive noise 
spikes and thus the measurement of the average submillimeter 
flux would be biased upward. To investigate the effect of the 
beam on the mean measurement we have performed this anal- 
ysis for the 850/im-detected radio sources in the CUDSS 14 
field (i.e. we measure the average submillimeter flux at the 
radio-determined position). The true average submillimeter 
flux of these systems is 5.4 ± 0.4 mJy and our recovered av- 
erage is 5.2 ± 0.4 mJy. Thus this technique recovers the true 
average flux within the uncertainties. 

There are two further complications in this analysis, that we 
have not yet considered. (1) Are the results strongly affected 
by confusion? (2) To what degree are these measurements 
dominated by a small number of very bright submillimeter 
sources? We have already addressed the issue of confusion 
between the EROs and those submillimeter-sources which are 
not EROs by removing their submillimeter flux from the map 
but we have not discussed confusion between EROs them- 
selves. If two or more EROs are separated by a distance 
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Table 1 . EROs with possible SCUBA detections at > 3cr 



SCUBA ID 


K 


(I-K) 


(R-K) 




Note 


14.1 


19.3 ±0.1 


5.1 ±0.3 


>6.3 


8.7 ± 1.0 


1.4GHz det. 


14.2 


19.6 ±0.1 


4.9 ± 0.4 


5.7 ±0.4 


5.5 ± 0.9 




14.5 


20.7 ±0.1 


4.1 ±0.5 


>4.9 


4.6 ± 1.0. 




14.7 


20.6 ±0.1 


4.1 ± 0.5 


>4.9 


3.2 ± 0.9 




14.8 


19.7 ±0.1 


4.3 ± 0.3 


4.8 ± 0.2 


3.4 ± 0.9 




14.9 


19.2 ±0.1 


> 5.7 


>6.3 


4.3 ± 1.0 


1.4GHz det. 


14.23 


20.0 ± 0.1 


>4.7 


> 5.3 


2.8 ± 0.9 




3.4 


19.8 ±0.1 


4.0 ± 0.2 


>4.6 


8.0 ± 1.0 




3.7 


20.4 ± 0.1 


> 5.0 


> 5.1 


8.2 ± 1.5 


1.4GHz det. 


3.22 


19.6 ±0.1 


4.2 ± 0.2 


5.7 ± 0.5 


3.1 ± 1.0 





roughly less than the beam size (~ 15") they will be con- 
fused and their measured flux densities will be incorrect. In 
the simplest case where two EROs lie directly on top of one 
another, each emitting lmJy at 850/im, we would erroneously 
measure 2 mJy for each. As we will discuss in more detail in 
§7 a number of the EROs associated with discrete SMGs have 
apparent ERO companions (§4) and therefore the naive treat- 
ment outlined above will lead to an over-estimate of the mean 
submillimeter flux. 

We now employ a more sophisticated statistical analysis, 
accounting for this confused flux in the following way. Con- 
sider two confused EROs with true submillimeter fluxes of I\ 
and h and separated by a distance r which is less than the ap- 
proximate beam size. The measured fluxes f\ and /2 may be 
written as follows: 

f^h+he^ 1 ' 2 ^ (1) 

f2 = h + he ( - r2 ^ ) (2) 

As we know f\, /2, and r we may solve for I\ and h- These 
equations are easily expanded for any number of confused ob- 
jects though in these data groups of more than three EROs are 
rare and there are no groups which contain more than four ob- 
jects. Using this procedure will result in an improved estimate 
of the mean submillimeter flux for the ERO sample as a whole 
though the measured flux densities of individual EROs may be 
incorrect since in reality only one of the confused EROs may 
be responsible for the submillimeter emission. 

We present the average submillimeter flux for the different 
samples, obtained using this method, in Table [2] The beam 
correction is important for 14% (11%) of the EROs in the 
(I—K) > 4 ((R-K) > 5.3) sample, and the corrections to their 
individual flux measurements vary from ~ 10 - 50%. For ex- 
ample, source CUDSS03.7 (shown in Figure 6) contains three 
EROs within 4" of eachother, and within the beam of a dis- 
crete SMG. The flux measurements for these three EROs be- 
fore the beam correction are 7.7 mJy, 7.3 mJy, and 6.6 mJy, 
with a non-weighted average of 7.2 mJy. After the beam cor- 
rection has been applied these fluxes become 6.2 mJy, 3.8 mJy 
and 3.8 mJy respectively, with a non- weighted average of 4.6 
mJy. Looking at the entire sample this method has lead to a 
decrease in the measured mean flux of ~ 20% (Table|2j. 

In Fig. [5] we show the mean 850/im flux the K-selected 
galaxies in both fields as a function of their color, for both 
color definitions. Clearly, the mean submillimeter flux is a 
strong function of color, with redder objects producing more 
submillimeter emission, regardless of the color definition. We 
also plot the results of W02, denoted by the filled triangles 



and same trend is seen, though W02 measure systematically 
higher flux above (I-K) ^> 3.5. 

The question of whether these measurements are dominated 
by a small number of bright objects remains. W02 concluded 
that their signal did not arise from a minority of objects; after 
removing the five most positive and five most negative mea- 
surements from their analysis they saw little change in the 
measured value. If we perform a similar analysis on these data 
our detection also remains significant. This does not mean, 
however, that the signal is not dominated by a minority of 
sources. Firstly, out of the few hundred EROs in our sample, 
ten are the best identified counterparts of SMGs and simply 
removing five of these does not fully address the issue of all 
ten dominating the emission. A more important point, how- 
ever, which we have already touched on and which we discuss 
later in §7, is that these submillimeter-bright EROs are often 
found in pairs or triplets with other EROs. Thus, though ten 
ERO (defined by (I-K)> 4) have been identified as SMG 
counterparts, 17 (again for (I-K) > 4) are associated with 
discrete SMGs. By associated we mean that they are found 
within the beam of a discrete SMG. The mean measurement 
is dominated by these EROs, which themselves are a minor- 
ity of the ERO population and which are associated with an 
even smaller number of discrete SMGs. As we show in Table 
|2j significant flux is not detected from EROs which are not 
found within the confused beam of a bright, discrete SMG. 

5. THE CONTRIBUTION OF EROS TO THE FIR BACKGROUND 
LIGHT 

Accounting for all of the objects which produce the cosmic 
background light is crucial for our understanding of the cos- 
mic star formation history. W02 estimate that EROs, with 
(I-K)> 4.0 and K' = 21.6 produce - 50% of the inte- 
grated background at 850/im, and that half of this energy is 
emitted by objects with K'< 20. Taking our mean 850/im 
flux, and the surface density of EROs in our fields we esti- 
mate that EROs produce 7-1 1% of the background at 850/im 
where the range in values encorporates the two ERO selection 
techniqu es we have used a nd t he values for the b ackground 
given bvlPu get etafl ill 9961) and lFixsen et alJ d!998h . Clearly 
a marked disagreement exists between these results and W02. 
The fact that W02 reach one magnitude deeper in /^'cannot 
resolve this disagreement since they estimate that half of their 
measured flux is emitted from objects with K'< 20. 

The difference in these two estimates stem from a larger 
surface density of EROs in the W02 fields and a higher mean 
850/im flux per ERO, and we discuss each in turn. The W02 
data cover 1 1 .5 arcmin 2 over three cluster fields and at K' = 
21.6 they find a surface density of 8 arcmin" 2 which is signif- 
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Table 2. The SubmillimeterFlux of EROs 



Field 


ERO criteria 


(simple) mean 


with beam 


"bright" EROs 








Sssotim (mJy) 


correction 


removed" 


both fields 


(I- 


K) > 4.0 


0.59 ± 0.07 


0.40 ± 0.07 


0.10 ±0.07 (17) 


both fields 


(R- 


-K) > 5.3 


0.61 ± 0.09 


0.56 ± 0.09 


0.22 ±0.09 (11) 


CUDSS+03 


o- 


-K) > 4.0 


0.36 ±0.11 


0.31 ±0.11 


0.16 ±0.12 (5) 


CUDSS+03 


(K- 


-K) > 5.3 


0.12 ±0.16 


0.11 ±0.16 


-0.05 ±0.17 (2) 


CUDSS+14 


(/- 


-K) > 4.0 


0.74 ± 0.09 


0.45 ± 0.09 


0.06 ±0.10 (12) 


CUDSS+14 


(R 


-K) > 5.3 


0.82 ±0.11 


0.76 ±0.11 


0.34 ±0.12 (9) 



'the number in brackets denotes the number of EROs removed from the analysis 



icantly higher than seen over larger areas ( Smi th etalJ20 02). 
Though their source counts are dominated by a single field 
when this field is removed the surface density still remains 
high (though within acceptable range). If the surface density 
measured by ( S mith et aT]l2002ft is used this reduces their es- 
timate of the background contribution by a factor of 0.7. We 
have covered roughly ten times the area of W02 and therefore 
our surface density should be a better reflection of the true 
surface density of EROs. 

The second issue is the mean 850/im flux per individual 
ERO. The largest difference here may stem from the effect 
of the beam correction (§4.2) which does not appear to have 
been included in the W02 analysis. In the case of our data this 
correction reduced the measured flux by roughly 20%. At 
our deepest limits we have a source density of ~ 2 arcmin" 2 
and, assuming no clustering, this corresponds to an average 
separation of EROs of roughly 48", and the beam correction 
is important for ~12% of the EROs. At the surface density of 
W02 sample of 8 arcmin" 2 the average separation of sources 
of about half this value (24") and the beam correction will be 
important for 45% of the EROs. If the highest density field 
of W02 is removed then about 19% of the sources will be 
effected. Thus, even for the lowest surface density used by 
W02 the correction is likely ~ 30% (and more than 50% if 
the highest surface density is used). 

We can also make some very simple arguments against 
such a large ERO contribution to the background light. We 
claim that the mean 850/im flux of EROs is dominated 
by a minority of discrete SMGs with Ssso/im > 3mJy. In 
blank field submillimeter surveys such objects produce ~15- 
30% of the FIR backgro u nd at 8 50itm jBlain et alJ 119991: 
Bareer^jgwie^fc Sanders 1999; ICowie. Bareer. & Kneibl 
120021 ISmail et alJ l2002bl IWebb et alJ l2003bl) and 20% 
of these objects have possible ERO counterparts to our 
magnitude-limits (§4). Therefore, a contribution on order 
6% from these EROs is expected. The lack of a significant 
detection of EROs below the 3 mJy means that the total 
contribution to the background cannot be much higher than 
this. 

In Fig.|4]we plot the mean 850/.tm flux of EROs in our sam- 
ple as a function of their K-magnitude, and over-plot the mea- 
surements of W02. W02 reported a strong increase in submil- 
limeter flux to brighter magnitudes which is not reflected in 
our data. This is difficult to understand as it implies too many 
K' -bright counterparts to discrete SMGs. Specifically, the av- 
erage 850^m flux of K'~ 19 EROs in the W02 sample is ~ 
3mJy and the integrated source counts of EROs at this mag- 
nitude roughly equals the integrated source counts of discrete 
SMGs with 585o^tm > 3mJy. This implies that a significant 
fraction of the discrete SMG population should have ERO 



counterparts brighter than K'~ 19, and it is well established 
that this is not the case jW£bjg_g^lJ|20_0_3d IClements et alJ 
l2003tlSmail et alJ2002bHlvison et alJ2002l) . This holds even 
if we take the more realistic view that the average flux mea- 
surement reflects a wide range in actual ERO submillime- 
ter emission. Thus, While it is clear that the W02 sample 
contains a number of submillimeter-bright EROs at bright K- 
magnitudes it is not reasonable to extrapolate this flux to the 
entire ^'-bright population. 

6. THE NIR COLOR-COLOR DIAGRAM 

In the above analysis we have treated the entire sam- 
ple of EROs as a uniform population while the well-known 
reality is the population is heterogeneous. The presence 
of old elliptical galaxies which have no substantial dust 
will dilute any s ubmillime ter sig nal from the dusty galax- 
ies. Pozzetti & Mannucci (2000) argued that old passively 
evolving stellar populations may be statistically distinguished 
from young dusty starbursts based on their locations on a 
NIR color-color diagram. The reason for this is straight- 
forward. Between redshifts z ~ 1 -2 the RIJK filter bands span 
the 4000A break seen in elliptical galaxies resulting in (ex- 
tremely) red (I-K) and (R-K) colors while the (J-K) colors 
remain moderate. Dusty starbursts have a relatively smooth 
SED in this region producing red colors in all filter combina- 
tions. The success of this technique to isolate dusty EROs, at 
least in a statistical sense, has also been reported by W02. 

In Fig. |5] we show the NIR color-color diagrams for the 
ERO samples, using z (9000A) in place of J. Also shown 
are the locations or limits of the EROs identified with dis- 
crete submillimeter sources (S^o^m > 3mJy). Overlaid are 
the tracks of two template galaxies ranging from z = 1 - 2: a 
starburst with dust and an early-type. The dotted line shows 
the rough division line between areas inhabited by each SED. 

This rough division can be used to group the EROs into 
those with elliptical-like colors and those consistent with star- 
bursts provided the redshifts are restricted to z = 1 - 2. Beyond 
this redshift old early-types will begin to contaminate the star- 
burst side of the plot. Though the error-bars are substantial 
the submillimeter luminous objects clearly lie on the expected 
star-burst side of the diagram. Ideally one would now repeat 
the mean 850/im flux measurement on each separate group 
to see if the dusty starbursts as a population indeed produce 
more submillimeter emission than the early-types. In this case 
however, the original measurement on the heterogeneous pop- 
ulation was entirely dominated by the discrete submillimeter 
sources which now almost exclusively lie on the starburst side 
of the plot. Thus, a repeat of this analysis would increase the 
submillimeter signal (by statistically removing old elliptical- 
like EROs) but would not provide any new information since 
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FIG. 4. — The mean 850/im flux measured for EROs as a function of their 
K magnitude. Top: The (/— ST)-selected sample. The solid circles correspond 
to this work and the solid triangles correspond to Wehner. Barger. & Kneib 
(2002). Bottom: The (i?-AT)-selected sample. 



we would still only be detecting these bright systems. 

7. ERO COMPANIONS 

It is well establishe d that t he ERO population is str ongly 
clustered JDaddi et alJ l2000t IRoche et alJ 120021 iFirth et alJ 
120021 IRoche. Dunlop. & Almainil 120031) over a wide range 
of magnitudes and color selection criteria, and this is easily 
interpreted in the framework of hierarchical models: galax- 
ies which are old at z ~ 1 (as a significant fraction of the 
EROs are) must have formed in the highest density regions 
in the early universe, and these regions themselves are highly 
clustered. Since separating the old elliptical-like galaxies 
from the starbursts is not trivial the contribution of each of 
thes e groups to the clu stering signal difficult is to disentan- 
gle. iDaddi et"al] J2002) however, have shown that while old 
ellipticals are indeed highly clustered, the dusty starburst pop- 
ulation is not JPaddi et alJ (12002). K < 19, {R-K) > 5). 
IRoche. Dun lop. & Almaini (2003) have investigated the fre- 
quency of close pairs of EROs ((I775 — > 3.92, Ks < 22) 
and found no evidence that dusty star-forming EROs are more 
likely to exist in close pairs than old passive EROs, or even 



optically selected galaxies dPatton et alJl2002l) . Thus, there 
is very little evidence for strong spatial correlations between 
dusty starburst EROs. 

The small areas observed thus far with SCUBA have 
not allowed a reliable clustering measurement for the SMG 
popul ation ( IScott et alJl2002t IWebb et al.ll2003bt iBorvs et alJ 
2003). These systems may be the result of purely stochas- 
tic processes and therefore would not b e clustered, as seen 
for the local FIR-se lected population ( iMoore et alJ 119941 
iKim & Sandersll998l) . An alternative scenario is one in which 
the SMGs are strongly clustered, and represent sign-posts for 
high-density, active environments at high-redshift. This latter 
picture is generally supported by the large mass estimates for 
a small number of these systems (>10 n M^~) (iFraver et alJ 
[T999tlTvison etZl20 00: Oenz el et all2 00l lNeri et alJ2003J) . 
Thus, though dusty EROs as a group may not be highly clus- 
tered, perhaps there is a subset of systems which are found in 
high-density regions and it is these which are related to the 
bright SMG population. If so, this implies environmental ef- 
fects are driving the star-formation in these systems. 

There are ten SMGs in the CUDSS catalog which are best 
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identified with EROs (see §4.1). Of these, five have another 
ERO within 10", and three of these five have two additional 
EROs within this distance (Fig. |6j. To quantify the signif- 
icance of this we note that a randomly chosen ERO in our 
sample will have a companion within 10" ~30% of the time. 
As there is no spectroscopic redshift information for these 
groups of EROs, it is impossible to know if all of these are true 
physical pairs or simply projections, and to determine their 
true separations. Nor can photometric redshifts constrain the 
redshifts with enough accuracy since all of these systems are 
very faint and the photometric wavelength coverage is insuffi- 
cient for reliable redshift estimates. However, there is circum- 
stantial evidence which suggests that most of these pairs are 



real. Four out of the five groups consist of EROs with similar 
apparent K magnitudes and NIR colors (the exception being 
CUDSS14.5). If the groups were due to the chance align- 
ment of galaxies at different redshifts such uniformity would 
not be expected. This argument is clouded by the fact that 
one of the EROs in each of the CUDSS03.7 and CUDSS14.9 
groups have radio emission, while their neighbors do not, 
which could indicate different redshifts. 

The offsets between these EROs range from ~ 2" to 10" 
which correspond to projected distances of roughly 20-100 
/r'Kpc at redshifts of z ~ 1-2. At these distances the sys- 
tems are unlikely to be strongly interacting (for example, in 
the late stages of a merger), but could be bound systems ei- 
ther influencing each other, or experiencing the same global 
environmental influences. Due to the faintness of the objects 
and the quality of the K data the morphologies are inconclu- 
sive: though some objects show asymmetries, there are no 
clear indications of merger-like activity. 

In Fig. we show the mean 850/^m flux of EROs as a 
function of their nearest neighbor distance. That is we mea- 
sure the mean submillimeter flux of those EROs with at least 
one companion within a given distance, for varying distances. 
We find EROs with companions mark on average regions of 
strong submillimeter emission. Though these measurements 
are dominated by the very bright groups discussed above and 
shown in Fig [6] there are very few close pairs of EROs which 
do not show evidence of any submillimeter emission. 

This result can be interpreted in two ways. Firstly, the sub- 
millimeter emission is produced or dominated by only one 
ERO within these groups, and this object serves as a marker 
for an over-density of EROs. This view is supported by radio 
detection of one ERO in two of these systems (CUDSS14.9 
and CUDSS03.7) which likely marks the true location of the 
submillimeter emission. The isolation of the submillimeter 
emission to one component of a composit e system has been 
seen b y other authors (e.g. Ilvison et"afl J2001I) : iNeri et all 
(2003)) though generally on smaller angular scales than seen 
here. Secondly, all of the EROs confused within the SCUBA 
beam may be responsible for the submillimeter emission each 
at a relatively equal level. Though individually the objects 
could likely lie below the 3cr detection limit of these sub- 
millimeter maps, the combined submillimeter emission from 
all would result in a single observed submillimeter source. 
This type of effect was clearly seen during the Mont e Carlo 
noise a nd source detection analysis of our maps in Eales et al. 
(2000). These two scenarios have different implications for 
the strength of the submillimeter emission from individual ob- 
jects, but both imply some level of interaction is driving up the 
average star-formation rate in these systems. 

The spatial association between bright SMGs and regions 
of galaxy over-densities has been previously observed for 
a number of systems identified with high star formation 
rates or AGN activity, includ ing Lyct emitters, radio galax- 
ies, EROs, and x-ray source s Jlvison et alJ l2000: Sma il et alJ 
20031 lAlexander etafl l2003 ), and indeed between these sys- 
tems themselves regardless of submillimeter information. In 
hierarchical formation models such regions, characterized 
by enhanced numbers of extreme objects, may subsequently 
evolve into the most massive galaxy clusters in the present- 
day universe. 

8. THE NATURE OF THE SUBMILLIMETER-BRIGHT EROS 

Ultimately, we wish to understand why some EROs are so 
luminous in the submillimeter while others show very little 



8 





p. I 

... i .... i 



..." 



t -::r.V;i%r&'. .'■<■■ 

- ,-i i a,-, .,- -i , --, =,--:u-.- ,: , ,r. 




17 18 19 20 

K magnitude 

FIG. 8. — A comparison of the NIR properties of the (possibly) 
submillimeter-bright groups of EROs with the general population of EROs. 
Solid points denote EROs which are associated with a discrete SMGs and 
which are found in pairs or groups of size less than the SCUBA beam. 
Crosses correspond to the entire X'-selected catalog (though only the red ob- 
jects are shown). Some of these SMG-EROs may not be true submillimeter 
emitters but may be due to chance alignments. Though very tentative, many 
of these EROs are brighter than the background ERO population, for a given 
color. 



FIG. 6. — The AT images (20" x20") of the five bright submillimeter sources 
identified with EROs that have apparent companions. Each image is centered 
on the ERO selected as the identification (from statistical analysis) and EROs 
are marked by circles. North is up and East is left. 
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submillimeter emission, if at all. Obviously, the most domi- 
nant and trivial effect is whether an ERO is a dusty star-burst 
galaxy, rather than an old early-type, but even within the star- 
bursts there is clear variation. We have seen in the previous 
section that environment plays a role, such that EROs within 
small groups or pairs often show strong submillimeter emis- 
sion. 

The strong increase in the mean submillimeter flux of EROs 



with NIR color seen in Fig. [3] could be indicative of a rela- 
tionship between submillimeter flux and redshift within the 
ERO population, since a given starburst SED will become 
increasingly red when placed at higher redshifts. Thus, the 
EROs which show submillimeter emission may lie at predom- 
inantly higher redshifts than those that do not. The median 
redshift of the su bmillimeter-selected populat ion of SMGs is 
z ~2-3 dSmail et all2002blllvison et all2002tlChapman et alJ 
2003b), and so the lower redshift sources in this distribution 
may over-lap with the high-redshift tail of this ERO sample. 
Of course, redder NIR color may simply be an indication of 
an increased level of extinction within these EROs, such that 
redder ER Os are more sub millimeter-bright because they are 
dustier (Chapman et al. 2003a|). 

In Fig. [S] we show the color and magnitude of EROs as- 
sociated with the discrete submillimeter sources. Since we 
plot all of the EROs within the beam of a discrete source it 
is possible that not all of these systems are true submillime- 
ter emitters. The bulk of the objects with measured colors are 
among the brightest objects in the field at each given color. 
Again, if we treat NIR color as a very rough redshift indicator 
then the submillimeter-bright EROs are intrinsically more lu- 
minous than average ^-selected EROs a given redshift. Thus 
not only are these systems currently experiencing higher than 
average star-formation rates but they may also have substan- 
tial populations of evolved stars compared to the rest of the 
starbursting EROs (provided they lie at z < 2). The fact that 
at least some of these systems exist in over-dense regions may 
indicate that these EROs have experienced a continual series 
of mergers (even the ones that currently appear isolated) over 
an extended time-scale, building up significant evolved popu- 
lations through a series so strong bursts. 

9. CONCLUSIONS 
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We investigated the submillimeter properties of EROs in the 
CUDSS fields and found the following: 

1 . We detect significant submillimeter flux from EROs in 
the CUDSS fields of 0.40 ± 0.07 mJy or 0.56 ± 0.09 
mJy for (I-K)> 4.0 and (R-K) > 5.3 ERO selec- 
tion criteria, respectively. This implies EROs to a depth 
of K < 20.7 produce 7-11% of the extragalactic back- 
ground light at 850/im. However, the flux is dominated 
by a minority of systems which are very bright at sub- 
millimeter wavelengths and overlap with the discrete 
population of submillimeter-bright galaxies discovered 
with SCUBA. 

2. EROs with strong submillimeter flux lie on the expected 
starburst side of the NIR color-color diagram, as sug- 
gested bv lPozzetti & MannucclfeOOOl) . 

3. Pairs or small groups of EROs mark, on average, re- 
gions of strong submillimeter flux, indicating that star 
formation rates are enhanced within these high-density 
areas. This implies, some fraction of strong submil- 
limeter emitters trace over-dense regions in the early 
universe. We are unable to tell from these data if more 
than one object within these groups is responsible for 
the submillimeter emission. 
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